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For quantum fields on a curved spacetime with an Euclidean section, we derive a general ex- 
pression for the stress energy tensor two-point function in terms of the effective action. The renor- 
malized two-point function is given in terms of the second variation of the Mellin transform of the 
trace of the heat kernel for the quantum fields. For systems for which a spectral decomposition of 
the wave opearator is possible, we give an exact expression for this two-point function. Explicit 
examples of the variance to the mean ratio A' = ^ \pyi P ' > °f the vacuum energy density p of a 
massless scalar field are computed for the spatial topologies of R d x S 1 and S 3 , with results of 
\Q ! A'(R d x S 1 ) = (d+ l)(d + 2)/2, and A'(S 3 ) = 111 respectively. The large variance signifies the 

importance of quantum fluctuations and has important implications for the validity of semiclassical 
0^ ' gravity theories at sub-Planckian scales. The method presented here can facilitate the calculation of 

stress-energy fluctuations for quantum fields useful for the analysis of fluctuation effects and critical 
' phenomena in problems ranging from atom optics and mesoscopic physics to early universe and 

Q black hole physics. 

^ ! I. INTRODUCTION 

■ Many important physical processes involve vacuum fluctuations of quantum fields. Famous examples are Casimir 
C\| | effect |lj and Lamb shift. These effects shed light on some basic issues of quantum field theory, and guide the 
. beginning probes and queries into establishing a viable theory of quantum fields. In curved or topologically non- 
trivial spacetimes (see, e.g., ||), basic issues such as the ambiguity in the definition of vacuum states, particles, 
energy, and the regularization of the energy-momentum tensor occupied the central attention of investigators in the 
70's. Amongst the many formalisms developed, the zeta-function |^,^| and point-separation regularizations are 
particularly relevant to our problem. In the 80's, two directions were noteworthy in this field: the backreaction 
of particle creation on the dynamics of the universe and the fate of black hole collapse ||, and the study of 
symmetry breaking and critical phenomena via interacting field theory in curved spacetime (see, e.g., ||). These 



I . studies are needed for understanding the dynamics of the inflationary universe, such as spinodal decomposition and 
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nucleation, and from it, astrophysical consequences such as entropy generation and galaxy formation. Since particle 
creation stems from amplification of vacuum fluctuations, and galaxy formation have them as seeds, both directions 
share a common need to understand better the properties of vacuum fluctuations of quantum fields and how they 
affect the dynamics of spacetime in the early universe and state of the matter at the observable classical domains or 
epochs. In the 90's, the theoretical underpinning of this issue was taken up afresh in the work of Ford and coworkers 
0,0, who investigated amongst other problems the effect of vacuum fluctuations in a quantum field on the causal 
structure of a quantum field theory, and in the work of Hu and coworkers ]T^ , |l3|] , who introduced non-equilibrium 
statistical mechanical concepts and techniques for the study of noise and fluctuations associated with a quantum 
field, with applications to backreaction problems in semiclassical gravity |l4|,[l5| and structure formation problems in 
inflationary cosmology Jl(]|l7j ■ 

In this paper, we aim to set up the basic framework for the calculation of the fluctuations of energy momentum for 
quantum fields in curved spacetimes, focussing on the generalization of the zeta-function regularization method flS] . 
Elegant and powerful, its application is unfortunately limited to spacetimes which admit an Euclidean section. But 
this includes many cases of physical importance, such as the de Sitter universe (in the 5 4 coordinatization), Kaluza 
Klein theory (M 4 x B D ), and finite temperature theory (in the imaginary-time formulation i? 3 x S 1 ). Later we intend 
to connect it with the point-separation method, which, as one of us has anticipated, actually contains untapped useful 
information about the statistical mechanical (or kinetic theory) properties of systems of particles interacting with 
fields in a curved spacetime. The result obtained here is also useful for the calculation of fluctuations of Casimir 
energy relevant in atom-optics [^9|, mesoscopic physics [Q, and correlations of quantum fields |^l|,^2| relevant to 
black hole fluctuations ]2^-p5t and possible Planck scale phase transitions j26| . 

The goal of this paper is to develop a regularized expression for the stress energy tensor two-point function of a 
quantum field in a curved spacetime. This is done by first expressing the two-point function solely in terms of the 
quantum field effective action and its variations with respect to the metric. The effective action is derived by using 
the trace of the heat kernel that corresponds to the field action. Since the effective action can be expressed as a 
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function of the Mcllin transform of the heat kernel, the needed variations of the effective action can be evaluated 
by considering the variations of the heat kernel. The variations of the heat kernel can be related to variations of an 
operator acting on the heat kernel in the varied heat equation. 

In Section II we relate the stress energy two point function to the second metric variation of the effective action. 
In Section III we show how to express this second variation in terms of the generalized zeta function for the system. 
Section IV presents the specialization of these expressions to the Klein-Gordon scalar field. In particular, we derive 
the expression for the variance of the energy density, and discuss the need for regularization. In Sees. V and VI we 
explicitly compute the energy density variance for a flat spacetime with one periodic dimension and for the Einstein 
universe respectively. We give a brief discussion of our results in Sec. VII. 



II. STRESS ENERGY TWO-POINT FUNCTION 



We start by considering the generating functional (or partition function) of a scalar field cj> in the Euclidean section 
S of a spacetime manifold M, 
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and its functional derivatives with respect to the metric: 
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where \0in, out > are the vacuua defined at the in and out states. In terms of the effective action W — logZ, this 
becomes 
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The expectation value of the quantum stress energy tensor is given by 
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In analogue with this, we define the correlation function for the stress energy tensor as 

(0,out|T ah (x)T cd (y)|0,in) 

SS 



^gjxj Sg ab (x)- 



(2.3) 



(2.4) 



(T ab (x)T cd (y)) 



(0,out|0,in) 

4 1 , , SS 
JO, out 



\Jg{x)g{y) z 



ab (x) 5g cd {y) 



|0,in) 



VgJxJgiy) 



S 2 W 



+ 



sw sw 



S 2 S 



Sg ab (x)Sg cd (y) Sg ab {x) Sg cd (y) x Sg ab (x)Sg cd (y) 



J 



(2.5) 



For any local action the last term will not contribute to the final expression for the stress-energy correlation function. 
For such an action, this expression will depend on x,y as a delta function 5(x — y). Thus it need not be considered 
when computing the correlation function for x ^ y. The autocorrelation is understood as resulting from this by taking 
the y — > x or coincidence limit. Recognizing the second term in the last line as a product of expectation values of the 
stress-energy tensor, we can define the bitensor 
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III. SECOND VARIATION OF THE EFFECTIVE ACTION 



The classical action of a scalar field 4>{x) is 

Sty] = \j duftx)H4>{x), (3.1) 

where if is a second order elliptic operator. From the spectral decomposition of this operator, H = A n |n)(n|, 
where n denotes the collective indices of the spectrum, the effective action can be expressed as 

W = -ilndet(H//i) = -^Tr ln^ = -^ln^, (3.2) 

n 

where we assume the zero modes of H have been projected out and /i is a normalization constant with dimensions 
of mass squared. This expression for the effective action is only formal since H is not trace class. We regularize the 
effective action and the expressions derived from it via the zeta function method HQ . 
The zeta function for this system is defined as 

(„(s) = Tre- si » H ^ = it s J2K s - (3-3) 

n 

For 2s > dim A4, this sum is convergent. Then an analytic continuation in s is found such that it includes a 
neighborhood of s = 0. From this analytic continuation the regularized effective action is given as 
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For positive H the definition of the gamma function yields 



= k'(0). (3.4) 
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(h(s) = =£t / t s - x Tr U (t)dt, U (t) = e~ tH , (3.5) 
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i.e. the zeta function is giv en as the Mellin transform of the trace of the heat kernel Uo(t). We know Tr Uo(t) ~ t~ d / 2 



for t —>■ 0. Hence for (3.5) to be convergent, we have the condition s > d/2. This is most often the expression from 
which an analytic continuation is derived. In fact, one of the main points of the zeta function idea is that formal 
expressions such as Tr Uo(t) need be modified by the introduction of a factor f ' . Then once the analytic continuation 
is found, one takes v — 0. 

We now consider the effect of two small metric perturbations 5i and 82 on the effective action. (They are assumed 
to be independent of the order with which they act). The response of the effective action to these perturbations is 

g + h + S2] - 
(Wff(s)) 
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To evaluate this, we use the Schwinger perturbative expansion |2S| ]. For U(t) = e~ t ( H+Hl > , where H\ = 5\H <C H, 

t 2 r 1 

Tr U(t) = Tr U (t) - Or [HxU^t)] + - / d Ul Tr [i?if/ ((l - «i)*)iTiUb(«it)] + ■ ■ ■ (3.8) 

z Jo 

Using this, we can write the response of the zeta function to these perturbations as 

&*i6r(«) = tjftt / dtf+1 I dui ^ TY [{SiH)U ((l - ui)t)(S 2 H)U ( Ul t)] 
Zl ( s ) Jo Jo 

+Tr [(6 2 H)U ((1 - mW^Hpoimt)]} . (3.9) 
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As it stands, it is not finite. When the traces are taken involving Uq((1 — u\)t) and Uo(u\t), the divergences at 
(1 — ui)t,uit — > are present. At this point, we modify the above expression for the second variation of the zeta 
function by introducing the factor [ui(l — wi)i 2 ] . At the end of the calculation, once the analytic continuation is 
found, we set v — 0. We view the introduction of this factor as an extention of the zeta function method to the 
situation where the second variation is needed. The replacement of Uo(t) — > t u Uo(t) is the spirit of the usual zeta 
function method and reproduces the usual results when applied to the traditional problems, such as finding the first 
variation, which produces the expectation value of the quantum stress energy tensor. After the change of variables 



u = (1 — u\)t, v = u\t 
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the twice varied zeta function transforms to 
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Considering the first trace in the above expression, we find 
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we can now write the stress-energy correlation bitensor as 
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where the s, v — > limit is understood. 

For the rest of this paper, we will specialize to the cases of homogeneous spacetimes. This implies AT^ bcd (x,y) 
will only depend on r = x — y. We find it convenient to average over all x. We can do this while leaving the points 
separated by r. Also, the homogeneity will usually lead to a degeneracy of the eigenvalues. Thus the collective 
quantum number n can be split into principal and degenerate parts n — > n, m and the eigenvalues only depend on 
n: X nm — » A„. This allows the sum over the degenerate indices to be done before evaluation of the zeta function. 
Putting all these together the stress energy two point function becomes 
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where O = J M dx, the volume of the manifold. If the manifold is noncompact, it is understood to be the unit volume. 

IV. FORM FOR THE KLEIN-GORDON FIELD 

We now develop the general form for the second variation of the zeta function for the Klein-Gordon field. In the 
Lorentzian sector, we use the MTW signature convention (— 1, 1, . . .), and thus in the Euclidean sector, we have the 
signature (1,1,...). We assume the metric can be given the form 
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where hij is the metric for the spatial section. We denote the time and spatial variables by x = (r, x), the invariant 
spatial volume form by dx, and the spatial manifold by £ (thus M. — S 1 x £). 
The wave operator for the Klein-Gordon field is given by 

d 2 

H = -U + ^R + m 2 = -— - s A + £i? + m 2 . (4.2) 

Let u n (x) be the eigenfunctions of S A: s A« m (x) = — k 2 m„(x), where n denotes the (collective) quantum numbers for 
the spatial part of the spectrum. We assume the u„(x) are orthonormal. The Euclidean time is made periodic with 
a peroid of /3 = 1/T, where T can be interpreted as a temperature. The eigenfunctions are thus given by 

e~ ik ° T 27T7io 

4>k .n{ x ) = ' m u n {x), k = ——, n = 0,±1,±2, ... (4.3a) 
VP P 

and the eigenvalues by (see e.g., p7|) 

A feo ,n = k 2 + K 2 n + £R + m 2 . (4.3b) 
From our definition of the stress-energy tensor ( |3.13| ) we find 

= -2V a ^V b + 5ah (V c ^V c + ^V c V c 0) + 2&<f>R ab . (4.4) 

This differs from the usual definition of T a b[i/j,<j)], but this is to be expected. The eigenvalues used are themselves 
different, since we have the extra r-degree of freedom. 

In this paper, we wish to concentrate on computing the autocorrelation of the energy density 

Ap 2 (x) ee lim ((p(x)p(y)) - (p(x)) (p(y))) , p = T 00 . (4.5) 

This will give a measure for the magnitude of the fluctuations of the stress-energy vacuum expectation value. This 
choice yields 

p[<f>, *!>]=- (d T (f>d T iP - ^d 2 T 4>) + (Vi^VV + ^ S A0) , (4.6) 
where i is summed over the spatial degrees of freedom only. This becomes 

p[<fi n ,<fi*n>] = ~ (k 2 + k k' + K 2 n ) 4> n <K,' + W„) (V 4 ^) . (4.7) 



when we specialize to eigenmodes (13). When we consider \p[<j) n , </>*/] | 2 the contribution from the </>* <j> n > Vi^V 1 ^*/ 
and its conjugate will vanish when summed over. This will be shown case by case. This assumption provides 

| P [^,C<]I 2 = {ki + k k' Q + n 2 n ) (k' 2 + k k' + kI) |0„| 2 icf + |(v^n) (w;,)| 2 ■ (4.8) 

At this point, we wish to take the zero temperature limit, i.e. — ► oo whereby Y^ =-oo (pV^ 71 ") /f^, dkQ. We 
can now do these sums/integrals easily, since they amount to calculating the moments of gaussians. We find 
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and there is zero contribution for any odd power of /co or /cq. The normalization of <f> n and n / provide a factor of [3~ 2 
which cancels the powers of [3 introduced in approaching the (3 — > oo limit. Introducing the functions 
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and 



^nn'(u,v) = E nn ,(u,v) + ® nn , (u, v) + Q n > n (v,u) + H nn ,(u,v), 
we define the zeta function 
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From it, we derive an expression for the regularized autocorrelation of the energy density as 
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where we have used 1/T (s) ~ s + 7s 2 + 0(s 3 ) (7 is Euler's constant) for s ~ 0. We have assumed v > d/2 + 1, in 
which case both C*(0, f ) and d(y(s, v)/ds\ s= o are finite. Thus to find the regularized expression for Ap 2 , we need to 
find the analytic continuation of C*(0, v) in v such that it is finite at v = 0. 



V. FLUCTUATIONS FOR £ = R D x S 



As the first application, we calculate in this section the variance of the energy density for a massless (m = 0) 
minimally coupled (£ = 0) scalar field on a (d + 2)-dimensional flat (R = 0) spacetime which is periodic on one spatial 
dimension with period L. For this geometry, the spatial eigenfunctions are 
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with I = 2tt/L. Denoting by x = (xi, . . . ,Xd) the coordinates for the open dimensions and z for the one compact 
dimension, we have S A = J2j=i ®x + ^2 an< ^ hence k 2 n = k 2 + l 2 n 2 , (k 2 = |k| 2 ). From this we see we should take 
k and n as the principal indices and the angular degeneracy of k as the degenerate indices, i.e. — > J dfld-i, 

integration over the unit d—1 sphere. Where there is no confusion, we will denote the volume of the sphere as S l . 
To evaluate (4.10) we sum over the degenerate indices and perform the volume average 
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from which (4.10a) becomes 
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Also, since VjUj-nVu^, = (kk' COS7 + l 2 nn')uk n u* k i n i where COS7 = k • k', the momentum correlation term (4.10c) 
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Here we have used J dfld-i J d£l' d _ 1 cos-f = and (A7) from the appendix: J dftd-i J dQ,' d _ 1 cos 2 7 = (S d ~ 1 ) 2 /d. 
Also, when summed over n and n', Qknk'n' vanishes. 

We now turn to the principal index sums for k and k'. First consider the case when n = (or n' = 0), this leads 
us to evaluate 
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Here (a, b) = (3/2, 0) or (1/2, 2), and 2a + b — 3 = 0. We assume a boundary on the non-periodic dimensions at large 
distance and moved to infinity. This is effected by having J °° dk — > lim e ^o / e °° • Then (|T^) becomes 
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where we have used v > d/2 in evaluating the k— integration. Thus we find there is no contribution from the n = 
or n' = terms in 

Turning to the case where n and n' do not vanish, we consider the function 
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In the integrand, there are terms with either factors of k d 1 or similiarly for kl . Using 
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we have the rule that upon doing the k, k' integrations, we get an overall factor of (r (|) /2) 2 (uv) d ^ 2 and each factor 
of k 2 becomes d/2u, and k' 2 becomes d/2v. Applying this rule, ( 4.11 ) has been determined 
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This form of the zeta function now allows us to perform the needed analytic continuation. We make use of the 
relation 
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where Cr( s ) is the Riemann zeta function. Recalling (4.13), the variance of the energy density is 
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Since T (— (d — l)/2) = — ((d+ 1)/2)T (— (d + l)/2), the second and third terms in the above expression cancel, leaving 
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by way of the Riemann zeta function reflection formula: 
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This is our result for the zero temperature variance of the vacuum energy density for a massless minimally-coupled 
quantum scalar field on a d + 2-dimensional spacetime periodic in one spatial dimension. To get a measure of the 
fluctuations, we consider the dimcnsionlcss quantity 
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For the system at hand, the energy density is pi 



S^ 1 /(A + 2 



(P) = -^ITT^r r ^— j &(d + 2). (5.17) 

Thus we get one of the main results of this paper: 

A>(Z = R^Si)= {d+1 f + 2 \ (5.18) 

Kuo and Ford jllj computed the same measure of the fluctuations for the case of E = R 2 x S 1 via a different 
method. Their result of A' = 6 is obtained here for d = 2. It is interesting to note that the relative amount of 
fluctuation increases quadratically with the dimension of the spacetime. 



VI. FLUCTUATIONS FOR E = S 3 



As a second example, we calculate the fluctuations of energy density for a massless (m = 0) conformally coupled 
(£ = 1/6) scalar field on a 3-dimensional space of constant curvature S 3 with radius a. The spacetime M. is then the 
Einstein Universe. 

We start by writing the spatial metric as Q 

3 

ds 2 = 7ab aV = ^Z> a ) 2 , (6.1) 

o=l 

where the a a, s are the basis one- forms on the three sphere satisfying the structure relations 

Ao a = \e a bc a b a c . (6.2) 

Here e% c , components of the totally antisymmetric tensor, are the structure constants for the rotation group 50(3). 
The 7 {, = l 2 S a i, are constants of the space (the principal curvature radii in a static mixmaster universe |30(] ) and 
for the Einstein Universe X = S 3 we have h = 1% = I3 = a/2. The curvature scalar is R = 6/a and the volume is 
fl = 2ir 2 a 3 . 

Using the Euler angle parametrization, the basis forms are 

o" 1 = — simpdO + cos ip sin 6d(p, 

a 2 = cos ipd9 + sin ip sin 6d<j), (6-3) 

a 3 = dip + cos 9 d<fi, 



< < 7T, < (f> < 2?r, < V < 4tt. 

Taking e a as the invariant vectors dual to cr° and defining the angular momentum operators L a = ie a , the spatial 
Laplacian becomes 



S A = J2 C 2 (e Q ) 2 = ~ (Li + L 2 2 + L%) = (6.4) 

0=1 



For the harmonic functions on S 3 we can use the SO(3) representation (Wigner) functions Dj CM (9,ip,<p), where 
./ = 0, i, 1, 1, . . . are the principal quantum numbers and K,M = — J, — J+1,...,J — 1, J are the degenerate 
quantum numbers. We will also find it convenient to use the SO(4) representation functions, the hyperspherical 
harmonics with principal quantum number n = J/2. From L 2 D J KA[ = J(J + 1)D J KM we find 



, 4J(J+ 1) n(n + 2) 2 tr> _{n + l) 2 



X n = < + £R = ± ^- (6.5) 
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along with the spatial eigenmodes 



ujkm 



n + 1 
Q 



(6.6) 



for S 3 



To compute S n „/, we first use the sum rule J2m" d mm"( d 'm> m"Y = $M,M' to S et J2mk d mk( d mk)* = 2 J + 1 



Since 



J dx.\u nK MU„'K' 



M' 



angular 



(n + l)(n' + l) 



f d X J2\ D KM\ 2 J2\ D 'K'M' 
J KM KM 



(4.10a) is given by 



(n + l) 2 (n' + 1 ) 2 
87r^ 2 (ut))i 



(n + l) 2 K + l) 2 

n 



1 n(n + 2) n'(n' + 2) 2n(n + 2)n'(n' + 2) 



uv a^v 



(6.7) 



(6.8) 



We now turn to the momentum correlation term Il„ n / . To facilitate the evaluation of the spatial derivative terms, we 
use properties of the generators of the Lie algebra, which for SO (3), is just the quantum theory of angular momentum 
(see e.g. pl[). We have 



/ 2 3 

-) i V(n + l)(n' + 1) L a (D J KM )L a (D J K r M ,) 



a J Q 

We also recast the spatial volume measure: 

3 p2ir pit 



M 



dx. = — 

4 



JO 



2tt 

sin^W / di/j = — I dto, 



(6.9) 



(6.10) 



where we use the notation of |31|] for dfl. We assume the integrand is invariant under -0 — > tp + 2n. This condition is 
satisfied by D J KM . 
Defining 



J J' 

E E 

K,M=~J K',M' = -J' 



a=l 



the momentum term (4.10c) becomes 



n„ 



(n + l)(n' + l) 4 
87rr2 2 (ww)i a% 2 



Returning to (3.11), we can express it in terms of the angular momentum operators as 

Bjji = ^ / 2 L + D KM L - D K'M' + 2 L - D KM L + D K'M' + L 3 D KM L 3 D 

KMK'M' J ~ 



J' 

K'M' 



(6.11) 



(6.12) 



(6.13) 



where L± = L\ ± iLi are the raising and lowering operators. Introducing the convenient notation L a , a € (+, — , 0) 
where Lq replaces L3 and the symbols 



C% = ^J{J + 1)-K(K±1), C° K = K. 
the action of the angular momentum operators on the harmonic functions is neatly given by 

r t-)J _ ( _;\—a/-i—a r\J j /-)./* _ _,— as-ia /-)./ 

^a^KM — \ L ) U K u K-a,M 1 a KM — 1 ^ K u K-Va,M- 



(6.14) 
(6.15) 







Introducing 

J' 



<a/3;6e>= ]T f dnL a Di M LpDi* M LsD^ M ,L e D^r M ,, (6.16) 

If (,f 7 r.'/ Jl,/7 7/ •* 



K,M=-J K',M' = -J' 



we can express (p.ll) as 



B.jj> = < 00; 00 > +i (< ++; > + < :++> + < +-; +- > + < -+; -+ >) 

+- (< 0-; +0 > + < 0+; -0 > + < +0; 0- > + < -0; 0+ >) . (6.17) 



Using ( 6.15 ), we can write ( 6.16j ) as 



C K a CKC K S ,C e K , j d^lD'l^ +j3 M D'l i * +< , M ,D J K _ a M D J K T_ S j 

M TCI M> J 



< a i 8e>=(-lf + H a +^+ t ^ Ci?C p K C$C K , / dnD^* +0<M D J K r +eM ,D J K _ aM D J K r_ SiM , (6.18) 

KMK'Mi 



The integral of the four-fold products of Wigner function is given by 

dtin Ji n J ' 2 n j3 * n Ji * 



J KM 

This can readily be seen by two applications of the sum rule 



■h 


■h 




Js 


Ja 




Mi 


M 2 


M ) 


V K 3 


K 4 





J 3 J4 J 



«- 2 H (v + 1) 2 P. J K 2. 1 m P. 2 < 1 it. m ■ < 6 ' M > 



JKM 

and the orthogonality property 



D J K \ Ml D£ M2 = £ (2.7+ 1) £ £ £ * * ^ D- M (6.20) 



dflD J K * M D J K , M , = -^-Sjj,5 kk ,5mm>- (6.21) 



2J+ 1 



Utilizing this result, ( 6.16 ) has the form 



J J' 

<af3;Se>=(-lf +e i a+fi+s+e 8ir 2 ^ J2 C K aC K C K' C K> (6-22) 

K,M=~J K',M'=-J' 

V- (0T»M1\( J J' J"\( J J' J"\( J J' J"\(J J' J" \ 

n , >\K + [3 K' + e K" )\K-a K' — 5 K" J \ M M' M" ) \ M M' M" ) ' 



J",K",M 

Using the orthogonality property of the 3-j symbols 



E [MM 1 mi ) ( M 1>V m 2 ) = (2Jl + l) ^^"n, 



(6.23) 



we can reduce the above four- fold product of 3-j symbols to a two- fold product by doing the M, M', M" sums to find 
the final form that is most useful to us: 

< a/3; Se >= (-lf +i 8ir 2 ^(2 J" + 1) (6.24) 

J" 



x 

KK'K" 



Uk u k^k> ^k> I K + p K > + e K" )\K-a K 1 -5 K 



10 



The triangularity of the 3-j symbols implies the condition a + (3 + (5 + e = for < a/3; Se > not to vanish. 
With this relation, we can evaluate the terms we need for ( 6.17 ): 



E 



KK'K 



KK' 



, J J' J" 



< ++: > = 



< +-;+- > 



< 



K K' K" 

^K^K^K'^K' 



KK'K" 

E 

KK'K" 

E 

KK'K" 



J J' J" \ f J J' J" 

K + l K' - 1 K" I [ K - 1 K' + 1 



--, , - I J J J 

^k^K' \ K-l K'-l K" 



/"f+^+ [ J J J 

^k^k' \ K + l K' + l K" 



<00;00> = 8tt 2 ^(2J" + 1) 

J" 

87T 2 ^(2J" + 1) 

J" 

< ++; — > 

8tt 2 ^(2J" + 1) 

J" 

> =8tt 2 ^(2J" + 1) 
J" 

< 0-; +0 > = 8tt 2 J2(2J" + 1) 
J" 

= < +0;0- > 

< 0+; -0 > = 8tt 2 ^(2J" + 1) 

J" 

= < +0;0- > 
We expect Bjji to assume the form 

Bjj, = 87r 2 (aJ 3 + 6 J 2 + cJ + d){af 3 + bJ 12 + cf + d). 
Determining the coefficients by evaluating Bjji for four different pairs of J, J', we find 



E KK ' C k C K ' ( K-l K' K") [ K K'-l K" 



KK'K" 



E KK'C+C+K' ( K + 1 K , k")\K K' + l K" 

KK'K" v ' v 



Bjj- = 8tt 2 J(J + 1)( J + 2) J'( J' + 1)(J' + 2) = — n(n + l)(n + 2)n'(n' + l)(n' + 2) 



and 



along with 



n„ 



(n + l) 2 (n' + l) 2 2 



n(n + 2)n'(n' + 2), 



87rfi 2 (ww)2 a 4 

(n + l) 2 (n' + l) 2 [ 1 n(n + 2) n'(n' + 2) 4n(n + 2)n'(n' + 2) 



uv arv 



87rfi 2 (uw)3 

Combining these results, the needed zeta function is 



1 ^ POO />oo 

n.n — 1 



C*(0,J/) 



87rft 2 



dv(uvY 



x (1 + (n 2 - l)u + (n' 2 - l)v + 4(?i 2 - l)(n' 2 - l)uv) , 



Using relation (5.12), the analytic continuation takes the form 



<*M = ^^; r(„-i) C,(2,-3) 2 



(6.25a) 

(6.25b) 
(6.25c) 
(6.25d) 

(6.25e) 

(6.25f) 

(6.25g) 
(6.25h) 

(6.25i) 
(6.26) 

(6.27) 
(6.28) 

(6.29) 



(6.30) 



+2r L-^\r(v+ £j Cfl(2^ - 3) (Cr(2i/ - 3) - Cr{2v - 1)) + 4r (v + l\ (Ca(2^ - 3) - ( R (2v - I)) 2 | (6.31) 
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Setting v = 0, we get finally the variance of the energy density of a scalar field on the Einstein Universe: 



Ap 2 



37 



(6.32) 



768007r 4 a 8 ' 



Using the result Q for the energy density, 



1 



(6.33) 



P 



480vr 2 a 4 



we find the dimensionless measure of the fluctuations in the energy density (5.16) is given by 



A'(S = S 3 ) = 111. 



(6.34) 



Thus for the Einstein Universe, the fluctuations in the energy density are indeed quite large. Effects due to the 
fluctuations of the metric will become important before the size of the Universe approaches that of the Planck scale. 



In this paper we have shown how the correlation function for the quantum stress energy tensor is related to the 
second metric variation of the effective action. This parallels the definition of the expectation value of the quantum 
stress energy as the first metric variation of the effective action. A physically meaningful expectation value is derived 
from a regularized or renormalized effective action. Likewise, the correlation function is defined here in terms of the 
second variation of the regularized effective action. The correlation of the stress energy tensor computed for two 
distinct points is finite regardless of whether the effective action is regularized or not (excluding lightlike seperated 
points for a massless quantum field) . It is only when the autocorrelation is computed that the issue of regularization 
arises. Nonetheless, we choose to use the correlation function defined in terms of the regularized effective action. This 
is a consistent approach since then it is defined as the second variation of the same object for which the expectation 
value is the first variation. 

Since we have only considered geometries for which a Euclidean section exists, we can regularize the effective action 
via the zeta function method. For a quantum system its generalized zeta function is given as the Mellin transform of 
the trace of its heat kernel. The effective action is then given by the derivative of the zeta function. This constituted 
our starting point for computing the second variation. The key to the zeta function method is to control the divergence 
of the heat kernel present when the Schwinger proper time vanishes. This is done by introducing positive powers 
of the proper time and then performing an analytic continuation in powers of the proper time. At the end of the 
calculation, the variable of analytic continuation is set to zero. This is consistent in that the introduced power can 
be relaxed to zero at any point of the calculation to recover the initial formal expression. 

The second variation of the generalized zeta function is facilitated by the Schwinger perturbative expansion, which 
shows how the trace of the heat kernel responses. Once we have this response, we only need the trace of a pair of heat 
kernels. To make this resultant expression meaningful, the key idea of the zeta function is: for each of these traces 
we introduce a power of the proper time variable for that trace. The stress energy correlation function is expressed 
in terms of traces of the system's heat kernel, regularized via the generalized zeta function method. This is one of 
the main results of this paper. The zeta function method allows us to relax the introduced powers at any time and 
recover the formal expression for the correlation function. 

For geometries admitting a mode decomposition of the invariant operator, we display the correlation function 
explicitly in terms of these modes. The specialization to homogeneous geometries is considered and the simplification 
this entails is explored. 

Resent results [ p~T] , p~9| have suggested that quantum fluctuations of the energy density may be significant for systems 
with non-zero vacuum energy density. Our result confirms this assertion and goes beyond. In particular, the variance 
of the energy density for a massless scalar field is found. (The variance is the coincidence limit of the energy-energy 
correlation function.) This measure of the quantum fluctuations is calculable since we have developed the correlation 
function in terms of the zeta function regularized effective action. Our results are in excellent agreement with the 
results of for Minkowski spacetime with one compact dimension. We have extended this work to flat spacetimes of 
arbitrary dimension with one periodic dimension and found that the variance grows quadratically with the dimension 
of spacetime. This may have unexpected implications for Kaluza-Klein theory. We also found the fluctuations for 
the Einstein Universe, which turns out to be more than ten times larger than the energy density. This shows that 
quantum fluctuations will become important at energy scales below the Planck scale and supports the suggestion 
pq| that critical dynamics at such scales could reveal interesting new phenomena. Knowledge of the higher order 



VII. DISCUSSION 
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correlation functions of the quantum stress energy may be necessary to account for the full backreaction effect of 
these large fluctuations of the quantum fields on the dynamics of the geometry (see, e.g., |^2[| ) and for investigating 
the issue of the viability of semiclassical theories f [l2||l3j ] at the Planck scale. 

Acknowledgement NP thanks Professor Elizalde for correspondence on the zeta function method. Part of the work 
was done when BLH was on leave at the Institute for Advanced Study, Princeton. This work is supported in part by 
NSF grant PHY94-21849. 
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APPENDIX A: CALCULATION OF / DQ D J DQ' D cos 2 7 

In this appendix, we calculate (3 denned by 

B d = [ dSl d f rf0^cos 2 7 = /3Vol 2 (5 d ) 
Js d Js d 



(Al) 



where (Kid is the volume measure on the d sphere S d and cos 7 = \x-x' |. Here x and x' are unit vectors in i? d+1 (i.e. 
they are 'points' on S d ). If we parameterize S d with the Euler angles Oi, i = 1, . . . d with < 0\ < 2ir and < 6j < n, 
j ^ 1, we have 



dQ d = sin' 



d-i , 



Using 



we get 



sin 



fc-i 



Voi^) = / 



• • • sm0 2 d9d ■ ■ ■ d6\. 

r(Mi) 

27r (d+l)/2 

„ dn <* = r (d+i\ ■ 
s d r (— ) 



Let x' be a point on S d and define a(x') = J dfl d \x-x'\ 2 . We have a(x') is independent of x' and hence 

B d = ( dVL' d a{x') = Vo\(S d )a(x ), any x e S d . 
is" 1 

We take xq = (0, . . . , a; ,d+i = 1 )• Then |x-x'| 2 = cos 2 9 d = 1 — sin 2 0^ and 



a(f ) 



d£ld 
= Vol(5 d ) 

= Vol(5 d ) 
= Vol(5 d ) 



d£ld sin 2 0^ 



i /; sin d+1 Mgrf 

r(f) rm 



d+ 1 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



Thus we get 



d+1 



(A7) 
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